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Abstract 

The cross section for bottom quark production is calculated for 
the HERA-B experiment. We consider both the order cross section 
and the resummation of soft gluon corrections in all orders of QCD 
perturbation theory. 
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1 Introduction 

The calculation of production cross sections for heavy particles in QCD is 
made by invoking the factorization theorem |l| and expanding the contribu- 
tions to the amplitude in powers of the coupling constant a s (fi 2 ). Recent 
investigations have shown that near threshold there can be large logarithms 
in the perturbation expansion which have to be resummed to make more 
reliable theoretical predictions. The application of these ideas to fixed-target 
Drell-Yan production has been the subject of many papers over the past few 
years 0. The same ideas on resummation were applied to the calculation 
of the top-quark cross section and differential distributions at the Fermilab 
Tevatron in [3-5]. What is relevant in these reactions is the existence of a 
class of logarithms of the type (ln(l — z)) l /(l — z), where i is the order of 
the perturbation expansion, and where one must integrate over the variable 
z up to a limit z = 1. These terms are not actually singular at z — 1 due 
to the presence of terms in 5(1 — z). However the remainder can be quite 
large. In general one writes such terms as "plus" distributions, which are 
then convoluted with regular test functions (the parton densities). 

In this paper we examine the production of 6-quarks in a situation where 
the presence of these large logarithms is of importance, namely in a fixed- 
target experiment to be performed in the HERA ring at DESY. This actual 
experiment has the name HERA-B || [7| and involves colliding the circulat- 
ing proton beam against a stationary copper wire in the beam pipe. The 
nominal beam energy of the protons is 820 GeV, so that the square root of 
the center-of-mass (cm.) energy is \/S = 39.2 GeV. Taking the 6-quark mass 
as mi, = 4.75 GeV/ c 2 then the ratio of mj/ y/S ~ 1/8. If we choose the renor- 
malization scale in the running coupling constant as m;, then a s (m 2 ) pa 0.2 
so a s (ml) h\(\fS/mb) ~ 0.4, indicating that perturbation theory should be 
reliable. 

In perturbation theory with a hard scale we can use the standard expres- 
sion for the order-by order cross section in QCD, namely 

a(S,m 2 ) = f dxi f 2 dx 2 J2 fii^i, fi 2 )fj(x 2 , ^(Tij^s = x x x 2 S, m 2 , /i 2 ), 

S Sx 1 ij 

(1.1) 

where the fi(x,fi 2 ) are the parton densities at the mass factorization scale 
fi 2 and the are the partonic cross sections. The numerical results for the 
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hadronic cross sections depend on the parton densities, which involve the 
choice of /z 2 ; the choice of the running coupling constant, which involves the 
renormalization scale (also normally chosen to be fi 2 ); and the choice for the 
actual mass of the 6-quark. In lowest order (LO) or Born approximation the 
actual numbers for the cross section show a large sensitivity to these parame- 
ters. The next-to-leading order (NLO) results follow from the work of the two 
groups H and || [Kj. However, even the NLO results do not completely fix 



the cross section. There is still a sensitivity to our lack of knowledge of even 
higher terms in the QCD expansion, which can be demonstrated by varying 
the scale choice up and down by factors of two. In general it is impossible to 
make more precise predictions for heavy quark production given the absence 
of a calculation in next-to- next-to-leading order (NNLO). 

In the threshold region one can improve on the NLO results. In this region 
one finds that there are large logarithms of the type mentioned above which 
arise from the soft-plus-virtual (S+V) terms in the perturbation expansion. 
These logarithms can be resummed to all orders in perturbation theory. We 
will see later that the gluon-gluon channel is the dominant channel for the 
production of 6-quarks near threshold in a fixed-target pp experiment. In the 
next section we give some results for the parton-parton cross section. Section 
3 contains the analysis of the hadron-hadron cross section which is relevant 
for the HERA-B experiment. We give results in NLO and after resummation. 
Finally in Section 4 we give our conclusions. 



2 Results for parton-parton reactions 

The partonic processes that we examine are 

i(k 1 )+j(k 2 )^Q( Pl ) + Q{p 2 ), (2.1) 

where i, j = g, g or i, j = q,q and Q, Q are heavy quarks (c, b, t). The square 
of the parton-parton cm. energy is s = (k\ + k 2 ) 2 '■ 

We begin with heavy quark production in the qq channel. The Born cross 
section in this channel is given by 

a^(s,m 2 ) = fa 2 s ^)K q ,NC F - s (3 (l + *f\ , (2.2) 
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where C F = (iV 2 — l)/(2iV) is the Casimir invariant for the fundamental 
representation of SU(N), K qg = N~ 2 is a color average factor, m is the heavy 



quark mass, \x denotes the renormalization scale, and (3 — Jl — Am 2 / s. Also 
N = 3 for the SU(3) color group in QCD. The threshold behavior (s — > 4m 2 ) 
of this expression is given by 



(0) / 2 



7ia 2 s (fx 2 )K gg NC F -p. 



(2.3) 



Complete analytic results are not available for the NLO cross section as some 



integrals are too complicated to do by hand. However in |10[ analytic results 
are given for the soft-plus-virtual contributions to the cross section, and 
for the approximation to the cross section near threshold. Simple formulae 
which yield reasonable approximations to the exact 0(a^) results have been 
constructed in 
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The analysis of the contributions to the gluon-gluon channel in NLO is 
much more complicated. There are three Born diagrams and different color 
structures. The exact Born term in the gg channel is 



ATra 2 (fi 2 )K gg NC F - jCV 
4m 2 



+ C* 



5 + 5t)' 




(2.4) 



where K gg = (N 2 — 1)~ 2 is a color average factor and Ca = N is the Casimir 
invariant for the adjoint representation of SU(N). The threshold behavior 
(s — > 4m 2 ) of this expression is given by 



a gg, thres( S 5 m 



za 2 a (fi 2 )K gg -NC F [4C F - C A }(3. 



(2.5) 



Again, the complete NLO expression for the cross section in the gg channel 
is unavailable but analytic results are given for the S+ V terms in |J . These 



were used in [|TTJ to analyze the magnitude of the cross section near threshold. 

In H an approximation was given for the NLO soft-plus-virtual (S+V) 
contributions and the analogy with the Drell-Yan process was exploited to 
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resum them to all orders of perturbation theory. We are discussing partonic 
reactions of the type i(ki) + j (A: 2 ) — > Q(pi) + Q(j>2) + g(k 3 ), and we introduce 
the kinematic variables ti = (k 2 — p 2 ) 2 — m> 2 , Ui = (hi — p 2 ) 2 — m 2 , and 
s 4 = s + ti + u\. The variable s 4 depends on the four- momentum of the 
extra partons emitted in the reaction and is zero for elastic scattering. The 
first-order S+V result in the MS scheme is 



2 d 2 a'-j\s,t l ,u 1 ) 
dtidui 



<7ij(s,*l,«ly 



X 



2Ca 



7Y 



Si 



A 



A m 2 
In — In — 



W - A) 



(2.6) 



where A is a small parameter used to distinguish between the soft (s 4 < A) 
and the hard (s 4 > A) regions in phase space. Here we define C q q = Cp, 
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Ca, 



<t^(s,*i,ui) = 7ra 2 s (n 2 )K qq NC F 



t 2 



+ 



2m 2 



(2.7) 



and 



<r M (s,fi,ui; 



2na 2 s (fi 2 )K gg NC F C F - C A 



x 



+ — + (l - 



m 2 s\ 



(2.8) 



We can define an analogous result for the qq channel in the DIS scheme (but 
not for the gg channel) which is 



d 2 a [l) 

'dM^ (M1 ' ul) = 

" 1 



2C f 



7Y 



a s (/i 2 ) 



x 



s 4 



In -4 + In — - A 



1 . 2 A , A , m J \ . . 



(2.9) 
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The resummation of the leading S+V terms has been given in B. The 
result is 



2 c/ 2 (T^ es (s,ti,Mi) 



dt\du\ 



afAs^Ux) 



dsi 
A m 2 

+ /(— — )*(**) 
mr /j, z 



6{s 4 - A) 



(2-10) 



where 



/ s 4 


m 2 \ 








In 2 4" 






= exp 




\m 2 


77 




7T 




m . 



[g4/wr] 

r(i + 7j) 



exp^-^j 



Expressions for A, a s , rj, and 7^ are given in B. 



(2.11) 



3 Results for bottom quark production at 
HERA-B 

In this section we discuss 6-quark production at HERA-B, and we examine 
the effects of the resummation procedure discussed in the previous section. 
Following the notation in || the total hadron-hadron cross section in order 
a k s is 

a%\S,m 2 ) = Y / I' dr^(T^ 2 )^\TS,m 2 ^ 2 ), (3.1) 

where S is the square of the hadron-hadron cm. energy and i,j run over q, q 
and g. The parton flux ^-(t, /i 2 ) is defined via 

^•(r,/! 2 ) = f 1 —Hijix, -,/x 2 ) , (3.2) 

Jt X X 

and Hij is a product of the scale-dependent parton distribution functions 
f^(x,fi 2 ), where h stands for the hadron which is the source of the parton i 

H ij {x 1 ,x 2 ,^)^ft l {x 1 ,ii 2 )f^(x 2 ,f J , 2 ). (3.3) 

The mass factorization scale /1 is chosen to be identical with the renormal- 
ization scale in the running coupling constant. 
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In the case of the all-order resummed expression the lower boundary in 
(3.1) has to be modified according to the condition so < s — 2ms 1 / 2 , where 
s is defined below (see 0). Resumming the soft gluon contributions to all 
orders we obtain 

a^(S,m 2 ) = ]T f 1 dr^^^a^rS^ 2 ,^), (3.4) 
ij Jto 

where is given in (3.24) of 0] and 

[m+ {m 2 + sq) 1 / 2 ] 2 
T o — g , w-oj 

with s = m 2 (/ig/ )i 2 ) 3 / 2 (MS scheme) or s = m 2 (/io//i 2 ) (DIS scheme). Here 
fjLo is the non-perturbative parameter used in H. It is used to cut off the 
resummation since the resummed corrections diverge for small /io- 

We now specialize to bottom quark production at HERA-B where y/~S = 
39.2 GeV. In the presentation of our results for the exact, approximate, and 
resummed hadronic cross sections we use the MRSD_' parametrization for 



the parton distributions [12]. Note that the hadronic results only involve 
partonic distribution functions at moderate and large x, where there is little 
difference between the various sets of parton densities. We have used the 



MRSD_' set 34 as given in PDFLIB [13[] in the DIS scheme with the number 



of active light flavors rif = 4 and the QCD scale A5 = 0.1559 GeV. We have 
used the two-loop corrected running coupling constant as given by PDFLIB. 
Note that we have checked scheme differences in the qq channel by using 
the MRSD_' set 31 as given in PDFLIB |L3| in the MS scheme. For the gg 
channel there is no DIS scheme so we always use the MS scheme. 

First, we discuss the NLO contributions to bottom quark production at 
HERA-B using the results in [8-10]. Except when explicitly stated otherwise 
we will take the factorization scale \x = where m& is the 6-quark mass. 
Also, throughout the rest of this paper, we will use m and mj, interchangeably. 
In fig. 1 we show the relative contributions of the qq channel in the DIS 
scheme and the gg channel as a function of the bottom quark mass. We see 
that the gg contribution is the dominant one, lying between 70% and 80% of 
the total NLO cross section for the range of bottom mass values given. The 
qq contribution is smaller and makes up most of the remaining cross section. 
The relative contributions of the gq and the gq channels in the DIS scheme 
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are negative and very small. The situation here is the reverse of what is 
known about top quark production at the Fermilab Tevatron where qq is the 
dominant channel with gg making up the remainder of the cross section, and 
gq and gq making an even smaller relative contribution than is the case for 
bottom quark production at HERA-B. The reason for this difference between 
top quark and bottom quark production is that the Tevatron is a pp collider 
while HERA-B is a fixed-target pp experiment. Thus, the parton densities 
involved are different and since sea quark densities are much smaller than 
valence quark densities, the qq contribution to the hadronic cross section 
diminishes for a fixed-target pp experiment relative to a pp collider for the 
same partonic cross section. 

In fig. 2 we show the K factors for the qq and gg channels and for 
their sum as a function of bottom quark mass. The K factor is defined by 
K = (a^ + er 1 ) | exact V ^ " 1 ? where is the Born term and | exact is the 
exact first order correction. We notice that all K factors are large. This is to 
be expected due to the new dynamical mechanisms which arise in NLO. The 
figure shows that higher order effects are more important for the gg channel 
than for qq. The K factor for the sum of the two channels is also quite large. 
However, the K factor for the total is slightly lower since we also include the 
negative contributions of the qg and qg channels. 

These large corrections come predominantly from the threshold region for 
bottom quark production where it has been shown that initial state gluon 
bremsstrahlung (ISGB) is responsible for the large cross section at NLO [|ll|]. 
This can easily be seen in fig. 3 where the Born term and the 0{a z s ) cross 
section are plotted as functions of r] cut for the qq and gg channels, where 
7] = (s — 4m 2 ) /4m 2 is the variable into which we have incorporated a cut 
in our programs for the cross sections. The cross sections rise sharply for 
increasing values of r] cut between 0.1 and 1 and they reach a plateau at higher 
values of r/ cut . This indicates that the threshold region is very important and 
that the region where s » 4m 2 only makes a small contribution to the cross 
sections. Note that in the last figure as well as throughout the rest of this 
paper we are assuming that the bottom quark mass is m\> = 4.75 GeV/c 2 . 

Next, we discuss the scale dependence of our NLO results. In fig. 4 we 
show the total 0(al) cross section as a function of the factorization scale 
for the qq and gg channels. We see that as the scale decreases, the NLO 
cross sections peak at a scale close to half the mass of the bottom quark 
and then decrease for smaller values of the scale. For the qq channel the 
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NLO cross section is relatively flat. The situation is much worse for the gg 
channel, however, since the peak is very sharp and the scale dependence is 
much greater. Since the gg channel dominates, this large scale dependence is 
also reflected in the total cross section. Thus the variation in the NLO cross 
section for scales between m/2 and 2m is large. 

In fig. 5 we examine the /io dependence of the resummed cross sections 
for 6-quark production at HERA-B for the gg and qq channels (in the MS 
and DIS schemes). We also show, for comparison, the /i dependence of 
pi®) + pi 1 ) | app -i-o-C 2 ) | app) where | app and | app denote the approximate 
first and second order corrections, respectively, where only soft gluon contri- 
butions are taken into account via the expansion of (2.11). Note that we have 
imposed the same cut on the phase space of S4 (54 > so) as for the resummed 
cross section. The effect of the resummation shows in the difference between 
the two curves for each channel. At small /io, a res diverges signalling the 
divergence of the running coupling constant. This is not physical and should 
be cancelled by an unknown non-perturbative term. There is a region for 
each channel where the higher-order terms are numerically important. At 
large values of /io the two lines for each channel are practically the same. For 
the qq channel in the DIS scheme the resummation is successful in the sense 
that there is a relatively large region of /io where resummation is well be- 
haved before we encounter the divergence. This region is reduced for the qq 
channel in the MS scheme reflecting the differences between (2.6) and (2.9). 
For the gg channel, however, this region is even smaller. 

From these curves we choose what we think are reasonable values for /i - 
We choose /i = 0.6 GeV for the qq channel in the DIS scheme (fio/m rs 13%) 
and /io = 1-7 GeV for the gg channel (fio/m « 36%). The values we chose 
for the qq and gg channels are such that the resummed cross sections are 
slightly larger than the sums + | app | app . Note that these /io 

values are not exactly the same as those used in ref. 4, where fio/m = 10% 
and fio/m = 25% for the qq and gg channels respectively, which predicted 
the mass dependence of the top quark cross section. In this reference the /i 
parameters were again chosen via the criterion that the higher order terms 
in the perturbation theory should not be too large. 

It is illuminating to compare fig. 5 with a corresponding plot for the 
top quark case (for instance with figs. 12, 13, and 14 in g, where the top 
quark mass was taken as 100 GeV). There one can infer that if we take the 
slightly larger /io values given above there is very little change in the top 
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quark cross section. The reason is that in this case the gg channel makes 
only a small contribution and the /io dependence in the qq channel reflects 
the small variation of the running coupling constant at a scale /i = 100 GeV. 
As the running coupling constant varies more rapidly at a scale /z = 4.75 GeV 
the Ho parameters should be taken from measurements at the lower scale and 
then used in the prediction of the top quark cross section. This emphasizes 
the importance of the proposed measurement at HERA-B. It is clear from 
fig. 5 that we cannot choose f^o/m = 25% for the gg channel for bottom 
quark production but we can choose Ho/m = 36% for the gg channel for top 
quark production, with very little change in the value of the top quark cross 
section. Both sets of parameters yield cross sections which are within the 
error bars of the recent CDF |T4]] and DO [JTJJ experimental results for the top 



quark cross section. Therefore our cut off parameters do have experimental 
justification. We would also like to point out that an application of the 
principal value resummation method has been recently completed by Berger 
and Contopanagos [|l6j leading to essentially the same mass dependence of 
the top cross section as reported in ref. 4, which again justifies our choice 
for /io- Finally note that we could just as easily have chosen to work in the 
MS scheme for both channels by changing /zo in the qq channel to /io ~ 1.3 
GeV. The reason the DIS scheme is preferred is simply because it has a larger 
radius of convergence. 

In fig. 6 we plot the NLO cross section for /z = m/2, m, and 2m, and the 
resummed cross section with the values of /zo that we chose from fig. 5, as 
a function of the beam momentum for 6-quark production at fixed-target pp 
experiments. The width of the NLO band reflects the large scale dependence 
that we discussed above. The total NLO cross section for 6-quark production 
at HERA-B (beam energy 820 GeV) is 28.8 nb for /z = m/2; 9.6 nb for \i — m\ 
and 4.2 nb for /z = 2m. The resummed cross section is 18 nb. This value 
was calculated with the cut s 4 > s while no such cut was imposed on the 
NLO result. As we know the exact 0(a^) result, we can make an even better 
estimate using the perturbation theory improved cross section defined by 

imp _ ^res , (1) I _n-W I ("? f,\ 

U H — U H ' U H lexact U H lapp ) l d - u / 

to exploit the fact that | cxact is known and a# | app is included in a r ff. 
Therefore we also plot the improved total cross section versus beam mo- 
mentum in fig. 6 (where we have taken into account the small negative 
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contributions of the qg and qg channels). The improved total cross section 
for 6-quark production at HERA-B is 19.4 nb. 

It is interesting to make an estimate of the theoretical uncertainty in 
the resummed cross section. For this one would have to determine the 
non-leading logarithmic factors in the order-by-order cross section. This is 
a difficult problem because the gg channel dominates so our knowledge of 
the order-by-order Drell-Yan cross section is not relevant. The problem is 
presently under investigation |T7|]. In [[| the logarithms involving the scale 



parameter fi were resummed into the factor exp(— 777^) in (2.11). This in- 
volves the assumption that these factors exponentiate, which may not be 
true. However we will use (2.11) as an estimate. The scale dependence en- 
ters into the running coupling constant, the parton densities, as well as the 
factors involving 77 in (2.11). To get some idea of the \x dependence of the 
resummed cross section in (2.11) we have rerun our programs with [i = 4m/ 5 
and n = 5m/4 and plotted the results in fig. 6. We see that the cross sec- 
tion is still sensitive to the choice of fi, mainly through the exponential in 
(2.11). As noted above this is probably misleading and the investigation in 
[Q should yield a more reliable result. 



4 Conclusions 

We have presented NLO and resummed results for the cross section for bot- 
tom quark production at HERA-B. It has been shown that the gg channel is 
dominant and that the threshold region gives the main contribution to the 
NLO cross section. The resummation of the S+V logarithms produces an 
enhancement of the NLO results and yields the value of 20 nb at = 39.2 
GeV for m;, = 4.75 GeV/c 2 . The theoretical error on this number is under 
investigation. 
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List of Figures 

Fig. 1. Fractional contributions of the gg (MS scheme, short-dashed 
line), qq (DIS scheme, long-dashed line), qg (DIS scheme, lower dotted line), 
and qg (DIS scheme, upper dotted line) channels to the total 0(af) 6-quark 
production cross section at HERA-B as a function of 5-quark mass. 

Fig. 2. The K factors as a function of 6-quark mass for 6-quark production 
at HERA-B for the gg channel(MS scheme, short-dashed line), the qq channel 
(DIS scheme, long-dashed line), the sum of the gg and qq channels (dotted 
line), and the sum of all channels (solid line). 

Fig. 3. Cross sections for 6-quark production at HERA-B versus r] cut with 
m b = 4.75 GeV/c 2 for the qq channel in the DIS scheme and the gg channel 
in the MS scheme. Plotted are the Born term (gg, upper solid line at high 
Vcut] QQ, lower solid line at high ?7 cut ) and the O(af) cross section (gg, upper 
dashed line; qq, lower dashed line). 

Fig. 4. The scale dependence of the NLO cross section for 6-quark pro- 
duction at HERA-B with m b = 4.75 GeV/c 2 for the qq channel in the DIS 
scheme (lower line) and for the gg channel in the MS scheme (upper line). 

Fig. 5. The /i dependence of the resummed cross section for 6-quark 
production at HERA-B with m& = 4.75 GeV/c 2 for the qq channel in the 
DIS scheme, in the MS scheme, and for the gg channel in the MS scheme. 
Plotted are a v q e ^ (lower solid line DIS scheme, middle solid line MS scheme) 
and a r g °g (upper solid line). Also we plot the sum + | app +cr( 2 ) | app 
(lower dotted line for qq in the DIS scheme, middle dotted line for qq in the 
MS scheme, and upper dotted line for gg). 

Fig. 6. Resummed, improved, and NLO cross sections versus beam mo- 
mentum for 6-quark production in fixed-target pp experiments for m& = 4.75 
GeV/c 2 . Plotted are the total resummed cross section (long-dashed line), the 
total improved cross section (short-dashed line), and the total 0(a^)) cross 
section (/i = m solid line, n — m/2 upper dotted line, \x = 2m lower dotted 
line). The scale variation in the resummed cross section with \i = 5m/ 4 is 
plotted as the upper dot-dashed line and that with \i = 4m/ '5 as the lower 
dot-dashed line. 
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